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TECHNICAL MEMORANDUM 


INSTRUMENTAL EFFECTS ON THE TEMPERATURE AND DENSITY DERIVED 
FROM THE LIGHT ION MASS SPECTROMETER 

INTRODUCTION 


In a recent article, Reasoner et ah [1] described the Light Ion Mass Spectrome- 
ter (LIMS) that was flown on the SCATHA satellite. This second- generation instru- 
ment was designed to measure the density, temperature, and flow velocity of low- 
energy ions (E < 100 eV). The instrument combined a retarding potential analyzer 
(RPA) with a magnetic mass spectrometer. Instrumental effects which are inherent 
in the design of the instrument are that the acceptance cone solid angle varies with 
mass and energy (Instrumental Effect A) and that the energy bandpass varies 
inversely as the mass (Instrumental Effect B). This report looks at how these effects 
influence flux measurements and derived temperatures and densities of the plasmas. 

A non-ramming, thermal plasma with a Maxwellian distribution is assumed in all the 
analysis. For the effects on the measured temperature and density, comparisons are 
made with the temperature and densities found using an instrument with a limited 
aperture (ft < 2ir steradian) but no acceptance cone energy dependence. 

The acceptance cone angle of the LIMS instrument decreases with increasing ion 
energy and mass. The effect is such that, generally, the density and temperature of 
the plasma are underestimated if a constant solid angle is assumed in the analysis. 
However, there is a small range near zero temperature in which neither the tempera- 
ture nor density is affected. 

The effect of the mass-dependent energy bandpass is to reduce the flux for the 
higher mass ions when their energies are close to the limit of the bandpass. 


PROBLEM DESCRIPTION AND SOLUTION 


The instrument is described in the article by Reasoner et al. [1]. The detailed 
description will not be repeated here, but it is noted that the LIMS is a planar-type, 
retarding potential analyzer (RPA) followed by a magnetic mass spectrometer. The 
field-of-view of the instrument is cone shaped with geometrically-set limits such that 
the half angle of the cone is a maximum of 60 deg. 

The energy passband of the instrument is determined by the fact that the mass 
resolution is given by: 

A M/M = 0.1 = C . (1) 

However, since E = kq/M, it follows that: 

A M/M = -A E/E = C , (2) 


Ci'" v V ‘ ** £ t - * 


./.:** ft ", ft 


and therefore 


i E = -CE = -CKq/M = C'/M . 


(3) 


nrtm rase is 

Ui ; POOR QUALITY 

Thus, as the mass gets larger AE decreases. If the bandpass is rectangular 
from 0 to the upper energy limit, E u , then for plasmas containing significant com- 
ponents with energy greater than E u , only part of the plasma will be detected by 
the instrument. If, in addition, there is a lower energy limit, E^, to detection se. 

by a retarding potential, as with an RPA device, then only that part of the plasma 
with energies between E^ and E y will be detected. This is shown schematically in 

Figure 1 in which the bandpass transmission, T(E), is drawn on a plot of the distri- 
bution function versus energy. The flux with energies between Ej. and E measured 
by the instrument can be expressed as: 


F (E l , E u ) = F (E l , co) - f (E u , «) , 

where the first argument of F, E L , is the lower energy limit and the second argu- 
ment, E y , is the upper energy limit of the flux integral, which will be discussed 

later. The fraction of the flux that is detected at an RPA equivalent to E r is 
therefore: 


R = F (E L , E u )/F (E l , “) = 1 - F (E u , «)/F (E l , ») . (4) 



Figure 1. Schematic illustration of Instrumental Effect B on the flux integral. 
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Equation (4) expresses the fraction of the flux that is measured in terms of an HP A 
expression for the flux into the instrument; i.e., F (E u , »). Note that equation (4) 

is not the fraction of the total possible flux but the fraction of the flux into the 
instrument. If there are no instrumental effects, these two should be equal. The 
results of the bandpass will be seen in a plot of flux versus RPA potential energy as 
a falloff in flux even when the plasma energy is less than the bandpass cutoff and 
the flux will approach 0 as E^ approaches E u . The effects of equation (4) on the 

RPA curves and therefore derived temperatures and densities will be discussed later. 

The behavior of the solid angle acceptance of the LIMS/SCATHA instrument is 
shown in Figure 2 from Reasoner et al. [1] (their Figure 7). From this gx'aph, it 
can be seen that a reasonable approximation to the behavior of the solid angle with 
inspect to plasma energy, expressed in terms of velocity, is an exponential of the 
form: 


n = n 0 exp (-v 2 ^ 2 ) , 


CBGSw; , , (5) 

OF POOR QUALITY 


where 

ft = solid angle of the acceptance cone for plasma energy E, 

ft Q = the solid angle determined by the geometrical limits of the instrument, 

V- = the characteristic velocity of the solid angle; i.e., the velocity at which 
the solid angle has decreased by a factor of e. 

varies with the mass as mentioned earlier. 



Figure 2. Computed values of the solid angle response ft in steradians as a 

function of ion energy and mass. 
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Assuming a non-drifting Maxwellian distribution, the flux into a planar-type RPA 
instrument is given in spherical coordinates by: 


G(V) 

Flux = (2nV t '/^if*) J J 


exp (-V 2 /V t 2 )(V 3 /V t 3 ) cosesineda d(V/V t ) , 


V=V 0 /oosG 


( 6 ) 


where 


Of KiOtf 


n = the density of the plasma, 


Vj. = the thermal velocity equivalent of the thermal energy of the plasma, 

V Q = the velocity required to overcome the retarding potential, 

0 = the angle between the normal to the instrument surface and the incoming 
velocity vector, 

6(V) = the maximum acceptance half angle for a given velocity, 

V = the velocity of an ion. 

The integration over the azimuth angle <f> has been carried out in equation (6) 
assuming azimuthal symmetry. 


A careful examination of the 0-V domain of the integral shows that there are 
regions of this space which are not accessible in the integral. Figure 3 shows these 
regions. Regions for which e > 9 max are not included in the integral due either to a 

limit imposed by the retarding potential or by the instrumental effect seen in the 

solid angle dependence on energy and expressed as 0 = G(V). Notice there is a 

max 

transition velocity, V', at which the limit on © max changes from that due to the 

retarding potential acting on the component of velocity parallel to the instrument 
normal to that due to the energy dependence of the solid angle. Reversing the order 
of integration, using the limits set as shown earlier and separating the integral over 
V into two parts, the flux integral becomes: 


Flux = (2n V /V7) 


V - e max= oos ( V V > 


/ / 

V=V 0=0 
o 


3 d(V/V t ) d0 


+ 


<» „max 

/ / [ 

v=y 0=0 


3 d(V/V t ) do 


(7) 
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where the integrand represented by the brackets remains the same as that in 
equation (6). 



Figure 3. The different regions of G - V space. Particles that have velocities 
with magnitudes and directions in the clear area contribute to the flux 
integral. Those that have velocities with magnitudes and direction 
in the stripped or stipled area do not. V 1 is the velocity 
magnitude at which the limits of integration change from 
being related to the RPA potential to being related 

to G(V) , 

Using the fact that 


8 = 2tt (1-cosG) or cosG ~ 1 - £2/2tt , 


with the solid angle of the acceptance cone determined by the particle's velocity, the 
expression for G(V) can then be written: 

G(V) = Q mnv = cos" 1 (1 ~ S2( V) / 2tt) 

III Cl A 

Using the variable change 


cose sine d0 = -cose d (cosO) 
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equation (7) can be written: 


OftIGS&A!. 

OF. P06R 


Qy^ilri 


V' 1 


Flux « ( 2n V 


t /* / 

*v J n 


oxp(-V 2 /V t 2 )(V 3 /V. t 3 )d(V/V t ) cosO d(cosQ) 


V J 0 =cos A <V /V) 
o max o 


+ (2nV 


t /*> /”/ 

V ' W^os’^HCV)/^) 


exp(-V 2 /V t 2 )(V 3 /V t 3 ) d(V/V t ) cosOd(cosU) . 


(8) 


Integrating over cos 6 gives: 


f V ' 

Flux = (nV t /^) I exp(-V 2 /V t 2 )(V 3 /V t 3 )(l-V 0 2 /V 2 )d(V/V t ) 


•V' 


+ (nV 


t //J) J expC-V S /7 t 2 )(V 3 /V t 3 )a-(l-Sl(V)/2 1 r) 2 )d(V/V t ) . (9) 


Whether the second part of this expression is analytically integrable or not depends 
on the form of fl(V). For the form chosen here [equation (5)], equation (9) is 
analytically integrable. 

Another change of variable is advantageous at this point. Let 
Z = V 2 /V t 2 = E/E t 

Z o = V o 2/V t 2 = V E t 

V = (V'/V t ) 2 = E'/E t 
K = nV t //7, = n/2E t /Mir , 


then equation (9) can be wri?+en: 


/* 2 ' KZ z" 2 ' r* 

Flux = T J exp(-Z)ZdZ + I exp(-Z)dZ+y/ exp(-Z)ZdZ + j[ exp(-Z) Z(l- JKZ)/2 tt) 

\ J v J v 


dZ 


( 10 ) 
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The transformation of Si(V) to 0(Z) is as follows: 


Qr, V fi'V * r-v 

yd'iUTY 


1 - 0(V)/2ir » 1 - (%/2 ir) exp (-V 2 ^ 2 ) = 1 - (fi 0 /2ir) exp (-Z V^/Vj*) 


2 nr 2> 


2 _ 


Let (V 1 /V t ) £ ‘ = Z x s E 1 /E t , then: 


1 - fi(V)/2ir si- (0 o /2tt) expC-Z/Zp = 1 - ft(Z)/2ir 


( 11 ) 


Using equation (11), expanding (1 - Si(Z)/2ir) , and integrating all tei’ms, equa- 
tion (10) reduces to: 


Flux 


_ K 


-( 


" Z o - Z « 

(Z o +l) e ° - (Z'+l) e L 


V \ -bZ 1 

•K-j (bZ’ + l) e 


Z o (e Z ° - e“ Z ') + (aZ’+l) e" aZ ' 


( 12 ) 


where a = 1 + 1/Z^ and b = 1 + 2/Z^. Notice there are six parameters on which the 
flux depends: E r JE’, E Q , E t , n, and M. 

In the limit of a constant Si (=S1 ), E- goes to infinity, a and b go to 1, and 
equation (12) reduces to: 


Flux ss 


K 


’ ” z n o -Z /cos 2 0 Z 1 

„ o _ 2 n o o 

e *- cos 0 e 
o 


(13) 


For limited but constant Si, E f can be found analytically to be 


E’ = E /cos 0 
o o 


and cos0 Q = (l~Si Q /27r) 


These last two equations were also used in reducing equation (12) to equation (13). 
Equation (13) can also be obtained by direct integration. For 0 Q = 90°, one obtains: 


Flux = | e~ Z o . 


(14) 


Since the LIMS is geometrically limited in the acceptance cone, all comparisons using 
equation (12) will be made against equation (13). 





INSTRUMENTAL EFFECT A 


OF PGUK QU&UI'Y 


From Figure 2, E^ for H + and He + ions was estimated to be 17 and 6 eV, 

respectively. These two ions were used so that comparisons could be made with the 
spacecraft-mounted instrument. Two examples of fits to data will be shown later. 
For the moment, however, this report will concentrate on the instrumental effect on 
derived temperatures and densities. 


In order to calculate the flux, it is clear from equation (12) that E r aa a func- 
tion of E q must be known. Because of the form assumed for ft, the equation for E' 

is a transcendental equation of the form: 


/E^7E T = l-(n 0 /2ir)exp(E'/E 1 ) . (15) 

The Z notation has been dropped in order to explicitly show the relationship between 
the energies. 

This equation is solved graphically and a table of values of E’ for the appro- 
priate values of E and Eh is generated. Figure 4 shows the graph used for H 
(E^ = 17 eV) to find the E r , one for each E Q . A graph such as this would have to 
be generated for each mass; since, as noted earlier, E^ varies with mass. Table 1 
gives the values of E q and E’ used in the calculations that follow. 



co, 0 


Figure 4. Graph of v^E /E q and l-( /2 tt) expC-F/E^) used to graphically find E f . 

For this set of curves E^ = 17 eV. 
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A set of curves for several different thermal energies, E t (=kT), was generated 

to see how the calculated flux varied as a result of the decrease in acceptance cone 
solid angle with energy (Instrumental Effect A). Figures 5 through 10 show the RPA 
curves, or plots of flux versus retarding potential, calculated for the LIMS using 
equation (12). Equation (lit) was used to calculate a reference set of curves which 
are plotted in Figures 5 through 10 as dashed lines. The solid lines are the RPA 
curves determined with equation (12) and the data in Table 1. RPA curves are 
+ 

shown for H (E- = 17 eV) in Figures 5 through 7, while Figures 8 through 10 show 

+ 

RPA curves for He (E 1 = 6 eV) . For each ion, three cases are presented: 


E t << El 


E f % E x 


and 


E. » Eh 
t 1 


cmamK & 

OF FOCI? QUALITY 


TABLE 1. VALUES OF E q AND E» USED IN THE CALCULATIONS 
OF FLUX IN EQUATIONS (12) AND (13) 



H + 

(E 1 = 17 eV) 

H e 

(E, = 6 eV) 

E 

E' 

E' 

0 



0.0 

0.00 

0.00 

0.2 

0.74 

0.36 

0.4 

1.40 

0.66 

0.6 

1.99 

1.16 

0.8 

2. oO 

1.59 

1. 0 

2.95 

1.98 

2.0 

5.01 

2.30 

4.0 

8.35 

3.73 

6.0 

11.00 

6.00 

8.0 

13.60 

8.00 

10. 0 

16.00 

9.81 

20. 0 

26.00 

11,60 

40. 0 

44. 00 

20.80 

60.0 

62.00 
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Figure 5. Plots of flux versus KPA for H + and = 1 eV. The solid line is the 

flux calculated including Instrumental Effect A while the dashed line is 
the flux calculated with the reference equation. 



e rpa 

Figure 6. Plots of flux versus RPA for H + and E^ = 20 eV. The solid line is the 

flux calculated including Instrumental Effect A while the dashed line is 
the flux calculated with the reference equation. 
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Figure 7* Plots of flux versus It PA for 11 * and ;; 50 oV . The solid lino is the 

flux calculated including Instrumental Effect A while the dashed lino is 
the flux calculated with the reference equation. 



*. 



Figure 8. Plots of flux versus BPA for Ho 1 * and E ( » 1 oV. The solid lino is the 

flux calculated including Instrumental Effect A while the dashed lino is 
the flux calculated with the reference equation, 


U 



e rpa 

+ 

Figure 9. Plots of flux versus RPA for He and E^. « 5 eV. The solid line is the 

flux calculated including Instrumental Effect A while the dashed line is 
the flux calculated with the reference equation. 



e rpa 

Figure 10. Plots of flux versus RPA for He + and E^ = 20 eV. The solid line is 

the flux calculated including Instrumental Effect A while the dashed line 
is the flux calculated with the reference equation. 
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Specifically, for H + (E^ = 17 eV) , E 

7, respectively; while, for He + (E ^ 
and 10, respectively. 


ORlQj)\f£r 

0F pooi? Ql^Xiry 

^ is set to 1, 20, and 50 eV in Figures 5, 6, and 
= 6 eV), E t is 1, 5, and 20 eV in Figures 8, 9, 


One fact immediately clear from the figures is that once the thermal energy of 
the plasma is greater than E 1 , the measured flux actually decreases with increasing 

thermal energy as a result of Instrumental Effect A. When the thermal energy is 
below E^, the flux seen by the instrument does not increase as rapidly with E t as 

does the flux affected only by a geometrically limited aperture. In fact, the flux 
increases relatively little in comparison to the reference. Also, the slope of the curve 
continuously changes with the retarding potential so that linear fits to the log of the 
flux are not appropriate and do not give valid results. 


A grasp of the magnitude of the influence of Instrumental Effect A on flux 
measurements can be obtained by comparing the flux calculated with a constant 
acceptance solid angle with the flux calculated taking into account the decrease in 
solid angle with energy. Equation (13) gives the flux for constant acceptance cone 
solid angle, while equation (12) gives the flux for a decreasing solid angle. Consider 
the ratio R of equation (13) to equation (12). This ratio for zero RPA gives the full 
effect of the solid angle decrease on the flux. The ratio at zero RPA is: 


R » (1 - cos 2 0 o )/(Q o a 2 /TT - ft Q 2 b 2 /4ir 2 ) 


(16) 


For the LIMS 6 =60°; therefore, ST = it, and equation (16) reduces to: 

O O 

R = 0„75/(a 2 ~b 2 /4) . 


(17) 


Table 2 gives the value of R for three ions for three values of E^. For the higher mass 

ions at the higher temperatures the reduction in flux due to Instrumental Effect A is 
quite dramatic. 

TABLE 2. RATIO OF FLUX AT ZERO RPA FOR NO EFFECT A TO FLUX AT 

ZERO RPA WITH EFFECT A 




H + 

Ile + 

0 + 

t 

(eV) 


E x = 17 eV 

E t = 6 eV 

> 

CD 

II 

T— 1 

w 


1 

1.08 

1.26 

1.42 



2.18 

5.89 

10.04 


20 

3.97 

15.31 
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INSTRUMENTAL EFFECT B 


0R1G&3 AL PA8F 2 Kj 

OF FOOK QUALIFY 

With the equation for the flux into the instrument [equation (12)1, return to equation (4) 
and calculate the effects of a limited bandpass (Instrumental Effect B). From the form assumed 
for the bandpass, an obvious effect of the bandpass will be the reduction of the integral flux for 
the energies between E^ and infinity, where the flux would fall to 0 for E L = E u . Figure 11 

shows the effect of bandpass widths of E u = 8 and 6 eV on a Maxwellian plasma with Ej. = 10 

eV. The acceptance field-of-view is assumed to be hemispherical. If this effect were dominant, 
it too could affect the shape of the RPA curve as shown and therefore the derived temperatures 
and densities. Since RPA devices detect integral flux, then for a plasma which would be 
affected (depending on the thermal energy and the bandpass), the RPA curve would not start at 
as high a flux as expected and would fall off faster than expected if the bandpass upper energy 
limit were greater than the plasma thermal energy. The bandpass effect for the Maxwellian 
plasma is obvious in Figure 11. However, the effects on an RPA curve may be much more 
subtle in instruments for which other instrumental effects are working; e.g., Instrumental 
Effect A. It should be obvious that for plasmas with E t « E u , there would be no effect on the 

RPA curve. Such effects due to the band-pass limitation would only be seen for plasmas that 
would have portions of the flux versus RPA curve above E u (E t ^ E u ) and in instruments for which 

E 1 1 E u* 


On the LIMS instrument, the bandpass has little effect on the flux for He and 

*|« 

H since E. << E for both ions and the solid angle decrease is the dominant effect. 
1 u + 

It would be significant for O for which E^ ^ 4 eV and = 6 eV. However, the 



Figure 11. Comparison of the RPA curve with and without bandpass limits. Two 
RPA curves with bandpass limits of E^ = 8 and 6 eV are shown. A Maxwellian 

plasma with E^ = 10 eV was assumed. No other instrumental 

effects are included in these curves. 
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flux of 0 into the instrument was never of such a magnitude that temperature or 
densities could be derived. Therefore, since that which follows is peculiar to the 
LIMS instrument, the bandpass effect will not be considered to be significant and is 
not included in the analysis. 


EFFECT ON PARAMETERS DERIVED USING CURVE FIT ROUTINES 


In a preliminary analysis of data from the LIMS, the function in equation (13) 
was used in a nonlinear least-squares fitting routine described by Bevington [2] to 
find the two parameters, temperature and density. To determine the error introduced 
by this procedure on the temperature and density derived, we used equation (13) to 
curve fit data calculated using equation (12). Such a procedure should show the 
effect of not including the decrease of the solid angle with energy (Instrumental 
Effect A) in the analysis. Figures 12 through 15 show the results of such fits as E, 
is varied. In these figures, the points calculated with equation (12) and used ns 
input data are shown as plus signs while the solid line is the curve calculated using 
equation (13) with the temperature and density to which the curve fit program con- 
verged in fitting the input data. In each figure, the input data is shown in the top 
line in the upper right-hand corner and the derived temperature and density in the 
bottom line. N is the input density, set to 10 in all the cases considered here, and 
M is the ion mass in amu’s. A somewhat extensive set of curves is shown so that 
the behavior of the RPA curve with E^ can clearly be seen. 

Both the thermal energy and density found with the curve fit routine rapidly 
diverge from the input values as the thermal energy increases. Figure 16 shows this 
result quite clearly. To generate this data, we used the derived temperature and 
density found as described earlier (Figs. 12 through 15). In the top panel of Figure 
16, the derived temperature is plotted as a function of the real or input temperature. 
For an input temperature of about 12 eV, the derived temperature would be about 
10 eV, The bottom panel of Figure 16 shows the ratio of the derived density to the 
input density as a function of input temperature. For example for an input tempera- 
ture of about 10 eV, the derived density is about half the real density. 

Good fits to the data can be obtained using equation (12) even with no correc- 
tion for drift velocity or spacecraft charge. Figures 17 and 18 show such fits for 
two different plasmas encountered on SCATHA. The solid line in Figure 18 is the 
best fit curve; the data is represented by circles. One of the reasons for the 
difference between the best fit curve and the data at the higher energy is simply 
due to the choice for the lower cutoff in counts. We have chosen it very near the 
noise limit and in fact are probably in the noise for counts less than about 10. 

Because of the effect on temperature and density shown in Figure 16, temperature 
and densities for the LIMS data are now derived by curve fitting equation (12) to 
the data. One must be cautious, however, in quoting temperatures and densities 
when the derived temperature of the plasma is larger than E^. When this condition 

is met, equation (12) becomes somewhat insensitive to temperature because now the 
solid angle decrease has an increasing effect compared to the RPA curve (compare 
the bottom two panels in Figure 15) . 
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RPA VOLTS RPA VOLTS 

Figure 14. Plots of fLux versus RPA potential for data calculated from equation (12) and fit with 
equation (13). The solid line is the best fit curve. The input and output parameters are shown 
in the upper right-hand corner of each panel. Here E (input) = 14, 16, 18, and 2i> eV. 
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Figure 16, Summary of the derived temperatures and densities versus E,. 
The top panel shows the derived temperature versus the input 
temperature. In the bottom panel is shown the fraction of 
the derived density to input density versus input 

temperature. 
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Figure 17. Example using equation (12) of a fit to data from the spacecraft 
instrument. This is an example of a cool plasma with - 0.76. The 

solid line is the best fit curve represented by E^. = 0.76 and N = 

11.56. The actual data is represented by the circles. 



Figure 18. Example using equation (12) of a fit to data from the spacecraft 
instrument. This is an example of a warm plasma with E^ = 8.34. The 

solid line is the best fit curve represented by E^ = 8.34 and N = 

3.124. The actual data is represented by the circles. 





SUMMARY 


An expression for the flux into a RPA is derived which takes into account the 
instrumental effect of a dependence on energy of the solid angle of the acceptance 
cone. A second instrumental effect of a limited bandpass is briefly discussed. It is 
shown that temperatures and densities derived without considering the effect of the 
solid angle dependence on energy will be too low, dramatically so for E t > E 1 , the e 

folding distance of the solid angle dependence. For E t << Ep there is effectively no 

impact on the derived temperatures and densities if the solid angle effect is ignored. 

Bandpass limits can alter the RPA curves further if both E t and E 1 are greater 
than the upper bandpass limit, E u . Thus, when either E 1 or E t approach E u> one 
must check to determine the effect of the bandpass limits. On LIMS, the bandpass 
would affect only the RPA curves for 0 + , since for 0 + ~ 4 eV and E u = 6 eV. 
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